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PACS 75 . 70 . Kw - Vortices in magnetic thin films 

PACS 75 . 78 . Fg - Magnetisation dynamics of domain structures 

PACS 75 . 10 . Hk - Classical spin models 

PACS 75 . 10 . Jm - Heisenberg model 

Abstract. - A two-dimensional easy-plane ferromagnetic substrate, interacting with a dipolar 
tip which is magnetised perpendicular with respect to the easy plane is studied numerically by 
solving the Landau-Lifshitz Gilbert equation. Due to the symmetry of the dipolar field of the 
tip, in addition to the collinear structure a magnetic vortex structure becomes stable. It is robust 
against excitations caused by the motion of the tip. We show that for high excitations the system 
may perform a transition between the two states. The influence of domain walls, which may also 
induce this transition, is examined. 
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Introduction. — Vortices in magnetic layers have 
been intensely studied [lJjT]- They are topological defects, 
and can be characterised by a chirality (an integer wind- 
ng number, which is a conserved quantity for topologi- 
cal reasons) and a polarity (the out-of plane magnetisa- 
tion). They can annihilate by antivortices having a wind- 
ing number of opposite sign. Vort ex- ant i vortex pairs are 
important excitations in two dimensional magnetic sys- 
tems. In open systems one can also excite isolated vor- 
tices, when the corresponding antivortex leaves the system 
at the boundary. 

The polarity of a vortex may be treated as a bit, as it 
possesses two very stable states which can be easily probed 
e.g. with GMR sensors, as those used in reading heads 
of magnetic disks. Thus magnetic vortex structures are 
promising candidates for novel non-volatile storage con- 
cepts. An important question in this context is, how to 
"write" a vortex. To switch the polarization, magnetic 
field pulses [8j[9], alternating magnetic fields 10 or spin- 
polarised currents [ll] have been proposed. 

In this Letter we present and analyse a new method, 
by which vortices can be generated or removed. The ex- 
citation energy is provided by a moving magnetic tip as 
used in magnetic force microscopy. We first show that 
vortex states are stable configurations in easy plane fer- 
romagnetic structures in the presence of a magnetic tip, 
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positioned above the ferromagnetic structure. Such a sys- 
tem has been studied recently in order to explore the fric- 
tion force decelerating the magnetic tip |12]jl4] . Then 
we show that the vortex remains stable when the tip is 
moved along the substrate, dragging the vortex through 
the substrate. Moreover the moving tip may also create or 
destroy a vortex structure, depending on the tip magneti- 
sation and velocity, offering an alternative way to switch 
a bit between effectively three states: A collinear state, as 
well as a vortex state with up or down polarity. Finally we 
study the stability of the dragged vortex structure when 
it passes through a domain wall. 

The system. — Our system consists of N=L x xL y 
classical Heisenberg spins Si=fjb i /jjL s on a square grid, 
where ji s is a material specific saturation magnetisation. 
The Hamiltonian contains two terms, corresponding to a 
substrate and a tip part 

T~L = ^sub + ^tip? (1) 

For the spin-spin interaction in the substrate we assume an 
isotropic exchange with interaction constant J and equiv- 
alent easy axes in x- and ^-direction, 

^sub = — J ^2 ' ^3 ~ d z ^2 ^i,z ~ ex V y^i^tx + $i,y)- 

i i 

(2) 

The easy-plane anisotropy d z <0 alone would lead to an 
infinite domain wall width, as all in-plane magnetisation 
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Fig. 1: (a) Sketch of the energy landscape. At least two local 
minima, corresponding to the CS (collinear state) and the VS 
(vortex state), exist. Which is the absolute minimum depends 
on the tip strength w. There is a potential barrier between 
the two states and thus the CS and the VS may be very stable 
against external driving (cf the blue arrows). A strong pertur- 
bation may lead to the transition from the VS to the CS (red, 
dashed arrow) and vice versa, (b-e) A part of a system near 
equilibrium, initialised in the CS (b,c) or the VS (d,e). The 
scanning velocities are v = 0.01 (b,d) and v = 0.2 (c,e). The 
colour coding as well as the arrows represent the magnetisation 
in the xy-plane, cf. fig. [4] 



configurations would be degenerated. Correspondingly, 
the Mermin- Wagner theorem would rule out long range 
magnetic order at finite temperatures 15 . The fourth- 
order anisotropy term e xy >0, however, breaks the contin- 
uous symmetry, so that the Hamiltonian Q has a fer- 
romagnetic low temperature phase with domain walls of 
finite width. We use the anisotropy parameters d z —— 0.1 J 
and e xy =0.1J in this Letter. 

The substrate-tip interaction is introduced via a dipolar 
term, 



3 (S^ • Gi)(Sti p ' Gi) — Si • Stip 



(3) 



where Ri = |R^| denotes the norm of the position of spin 



i relative to the tip 



r t i p , and e$ its unit vec- 



tor = Hi/Ri. Ti and r t i p are the position vectors of the 
substrate spins and the tip respectively, w is a free param- 
eter that quantifies the dipole-dipole-coupling between the 
substrate spins and the tip, thus controlling the strength 
of the tip. We use S t i p = (0,0, —1), w is a free parame- 
ter representing the tip strength. The tip is moved with 
constant velocity (v, 0, 0) two lattice constants above the 
substrate, along the middle line between two spin rows. 

Two different kinds of boundary conditions will be used. 
First we consider the case that all inhomogeneities in the 
system are due to the tip. Then it is natural to describe 
the system in the comoving frame of the tip 
y— direction, the boundaries are open, in the x 
dynamical. When the tip advances by exactly one lat- 
tice constant, the foremost row is duplicated, and the last 
one is deleted. That way, arbitrarily long times can be 
simulated. In the second part we investigate, how vortex 
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Fig. 2: Nonequilibrium steady state diagram with the pure 
states (collinear state, CS, and vortex state, VS), the co- 
existence state (CS & VS), the periodic state (PS) and the 
vortex-antivortex state (VAVS). The grey points represent the 
underlying simulations. 



generation is influenced by a domain wall that is pinned 
far away from the tip. In this case, a combination of 
open (x— direction) and fixed (y— direction) boundary con- 
ditions is more appropriate. The system sizes used are 
64 x 48 for the case of simulations in the absence of do- 
main walls, and 200 x 48 for the case where domain walls 
are of interest. 

The equation of motion is the Landau-Lifshitz-Gilbert 



(LLG) equation 16 
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dt 



[Si x hi + a Si x (Si x hi)} , (4) 



with saturation magnetisation /i s , gyromagnetic ratio 7, 
the phenomenological damping constant a (we use the 
high damping value a=0.5 in this Letter to reach a steady 
state in a short simulation time) and the local field = 
—dl-L/dSi. It produces Larmor precession with frequency 
I ^ I 7//x s , and a damping in the direction of the local field. 
In equilibrium, each spin points in the direction of its lo- 
cal field. To solve the LLG we use the Heun integration 
scheme [l8] . 

Nonequilibrium steady states. — Before studying 
the driven system, let us discuss the equilibrium case 
(v=0). If the tip is absent (w=0), the equilibrium con- 
figuration corresponds to all spins pointing in the same 
direction, the collinear state (CS). 

As we increase the tip strength, a second minimum ap- 
pears in the potential landscape, the cylindrically sym- 
metric vortex state (VS, cf fig. [T| (d)). At the same time, 
the CS minimum is continuously moved to the configu- 
ration of a slightly disturbed collinear configuration (cf 
fig. [l](b)). For any reasonably small tip strength, both 
equilibrium states, the CS and the VS, are stable. If we 
initialise a VS, and switch off the tip, the VS is still sta- 
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Fig. 3: Some snapshots of a system in the PS in a co- moving 
frame of reference, the tip positioned at the origin (the two lines 
correspond to the x- and 2/-axis). Vortex cores are marked as 
black dots, antivortex cores as white dots. Here we observe 
the periodical transition from the VS to the CS (a-c) via the 
release of the vortex from the tip, as well as the transition from 
the CS to the VS (c-f) via a VAVP creation. See the complete 
cycle as a movie online. 



ble: Being a topological defect, its winding number is a 
conserved quantity which cannot spontaneously change. 

Let us now move the tip with constant velocity (v>0): 
The system is driven into a nonequilibrium steady state 
(NESS), which is near the equilibrium for small excita- 
tions. Accordingly, in the NESS diagram (cf. fig. |2| we 
see at v>0 a region where the CS and the VS coexist. 
Here the initial configuration determines the finally sta- 
bilised NESS. As the excitation overcomes a threshold by 
increasing the tip strength w or scanning velocity v, the 
system may perform a transition from the VS to the CS 
(when the CS represents the total minimum) or vice versa, 
as sketched in fig. [TJa) by the red dashed arrow. In the 
NESS, we then observe either the pure CS or the pure VS. 

Another mixed state is the periodic state: Here the sys- 
tem is excited so strongly that neither the VS nor the CS is 
stable. The system then flips between the states back and 
forth continuously. We explain this state in more detail, 
because it sheds also more light onto the transitions to the 
pure states mentioned above. Let us start in the VS. Here 
the vortex is bound by the tip, because the cylindrically 
symmetric structure minimises the tip energy. However, 
in the PS (or when the CS is the corresponding NESS) the 
tip pumps so much energy into the substrate that the vor- 
tex may decouple from the tip, and moves away towards 
the system boundary (fig. [3jb)), leaving a CS behind (fig. 
[3jc). Accordingly the potential barrier in fig. [TJa) corre- 
sponds to the energy of a free moving vortex and a CS. If 
the system's NESS is the CS, the system has now reached 
its steady state. If the initial state was the CS, but the 
corresponding NESS is the VS or the PS, the system is 
again excited further by the tip, and a vortex- antivortex 
pair (VAVP) may nucleate under the tip (fig. [3^d)). The 
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Fig. 4: Legend for the colour coded plots, and the equilibrium 
configuration of a system with two domain walls, as it is used 
for the simulations as an initial state. The width of the middle 
domain W is a free parameter. 



energy of a VAVP represents the potential barrier which 
the system has to overcome at this transition. After the 
nucleation, the antivortex moves away from the tip. It 
may stay at a constant distance (e.g. fig. [3^e)), and we 
end up in a NESS called vortex-antivortex state (VAVS). 
If the antivortex moves out of the system, we again get a 
VS(fig.§f)). 

The two transitions lead to very different paths in con- 
figuration space. The reason why two different barriers 
have to be passed, as well as why one barrier is higher than 
the other, can be found in the topology of the system: For 
the transition CS— ^VS, a vortex must be created. How- 
ever, an isolated vortex cannot be created, as it represents 
a topological defect, which violates vorticity conservation 
(the total number of vortex and antivortex cores). Only 
VAVPs can be created, and thus the energy barrier here 
is higher than that of the VS— ^CS transition. The reason 
why the latter transition (which violates vorticity conser- 
vation) can occur at all, is that here the antivortex core 
interacts with the open system boundary. The different 
transition paths of the PS lead to a hysteresis. 

Influence of domain walls. — In the section above 
we studied a perfect system in order to get an insight into 
the occurring NESSs. However, a real magnetic specimen 
contains domain walls, which the tip has to pass through. 
The open question is, how stable the above characterised 
NESSs are against e.g. domain walls. In order to answer 
this question we now initialise a system which contains 
two 7r/2 domain walls in equilibrium, where the distance 
between the two domain walls (or the length of the middle 
domain) is a free parameter W (cf. fig. [4]). Additionally we 
fix the boundaries in ^-direction by the equilibrium mag- 
netisation, as it is sketched in fig. [4j to emulate infinitely 
long domain walls. The initial condition corresponds to 
the simulations initialised in the CS above, as we always 
start in the left domain. Let us start with a weak exci- 
tation (w—>0). Then the domain wall configuration is not 
disturbed by the scanning tip. We may observe a slight 
bending when the tip just passes by (like in the very left 
panel of fig. [5J, which soon relaxes again, like a rubber 
band. After the tip has passed the domain walls the con- 
figuration is again in its initial state. 
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Fig. 5: (As a movie and colour online) Fast magnetic tip (marked by the dashed circle, w = 3,v = 0.5), interacting with two 
domain walls, with the domain wall distance W = 20. The snapshots show a part of the whole system containing 50 2 spins. 
The high excitation may lead to the creation of a VAVP, see (c). As the NESS of the system is a CS, the VAVP is annihilated 
again, leaving the initial domain wall configuration. 




Fig. 6: (As a movie and colour online) Slow, but strong magnetic tip (w = 5, v = 0.3), interacting with two domain walls. Again 
only a part of the system (70 x 50 spins) is plotted. A VAVP is created. As the NESS for the present v — w combination is the 
VS, the vortex stays bound by the tip, and the antivortex travels out of the system, leaving a VS instead of the domain wall 
state. 



When the perturbation by the tip is stronger, a VAVP 
may nucleate at the domain walls, even if the energy of 
the moving tip alone is not sufficient, as the domain walls 
provide additional energy (cf. fig. [5J. Which state is fi- 
nally adopted depends on the NESS. If the corresponding 
state is the CS, the VAVP get annihilated, and the ini- 
tial domain wall state appears again. If the corresponding 
state is the VS, the distance between the vortex and the 
antivortex increases. Finally, the antivortex stays at the 
system boundary (as we have fixed boundary conditions 
in this part the antivortex cannot leave the system), leav- 
ing a VS that trails along with the tip. We observed that 
in the coexistence regime (CS & VS) the system ends up 
in the VS. The PS is not influenced by the domain wall. 
It may occur that at the domain walls additional VAVPs 
nucleate, which after a short lifetime annihilate again. 

The influence of the domain wall distance W is the fol- 
lowing: The smaller W, the larger is the energy density, 
the domain wall may provide to a VAVP creation process. 
Accordingly, the lowest tip strength w, at which a VAVP 
is created, is larger for larger W. 

Finally we discuss the influence of the boundary condi- 
tions. For instance from the right panel of fig. [6] one may 
claim that the fixed boundaries have a strong impact on 
the stability of the VS, and enforce the annihilation of the 
antivortex, as the separation of the VAVP generates con- 
tinuously growing domain walls. These may become arbi- 
trarily large and thus energetically much more unfavorable 
than the domain wall state. To get an annihilation pro- 
cess, vortex and antivortex must first move toward each 
other. This is only possible when the vortex is released 
from the tip first, as a free moving vortex or antivortex 
cannot move at velocities comparable to that of the tip. 
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Fig. 7: Influence of domain walls on the NESS. Empty circles 
represent simulations, where the system does not create any 
VAVP, filled circles represent these simulations where at least 
one VAVP is created. The colour coding denotes which final 
state has been observed: Blue circles represent the CS, yellow 
ones the VS, and green ones the PS. 

If the vortex is released or not from the tip depends solely 
on the current NESS, and not e.g. on the tailing domain 
walls, as they cannot provide additional energy density un- 
der the tip for a vortex release event. All observed VAVP 
annihilation processes occur directly after the tip-domain 
wall interaction, where the events under the tip may be 
seen as completely decoupled from the walls. 

Conclusion. — A tip, which is aligned above a fer- 
romagnet with an easy plane anisotropy, and magnetised 
perpendicular to the easy plane, energetically stabilises a 
vortex state in the substrate. The vortex state is stable 
against slight perturbations, occurring e.g. when the tip is 
moved with constant velocity v parallel to the substrate. 
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At a threshold velocity, which depends on the tip mag- [18] Garcia- Palacios J. and Lazaro F., Phys. Rev. B, 58 
netisation, the substrate may perform a transition from (1998) 14937. 

the vortex phase to the collinear phase and vice versa, 
but it may also exhibit a periodic switching for strong ex- 
citations. In summary we have three possible states the 
system may adopt (the CS state and the VS state with up 
or down polarity), and which may be switched by a mov- 
ing tip. The presence of a domain wall effectively shifts 
the state separation line between the coexistence state and 
the pure vortex state to lower tip magnetisation values. 
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